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Abstract. The galactic superluminal microquasar GRO J1655-40 started a new outburst in February 2005, after 
seven years in quiescence, rising to a high/soft state in March 2005. In this paper we study the X-ray spectra 
during this rise. We observed GRO J1655-40 with XMM-Newton, on 27 February 2005, in the low/hard state, and 
on three consecutive days in March 2005, during the rise of the source to its high/soft state. The EPIC-pn camera 
was used in the fast-read Burst mode to avoid photon pile-up. First, we contributed to the improvement of the 
calibration of the EPIC-pn, since the high flux received from the source required some refinements in the correction 
of the Charge Transfer Efficiency of the camera. Second, we find that the X-ray spectrum of GRO J 1655-40 is 
dominated in the high/soft state by the thermal emission from the accretion disk, with an inner radius of 13- 
14[D/3.2 kpc] km and a maximum temperature of 1.3 keV. Two absorption lines are detected in the EPIC-pn 
spectra, at 6.7-6.8 and 7.8-8.0 keV, which can be identified either as blended Fe XXV and Fe XXVI Ka and K/3 
lines, or as blueshifted Fe XXV. We find no orbital dependence on the X-ray properties, which provides an upper 
limit for the inclination of the system of 73°. The RGS spectrometers reveal interstellar absorption features at 
17.2 A, 17.5 A (Fe L edges) and 23.54 A (OI Ka). Finally, while checking the interstellar origin of the OI line, 
we find a general correlation of the OI Ka line equivalent width with the hydrogen column density using several 
sources available in the literature. 
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1. Introduction 

Galactic microquasars are accreting binary systems eject- 
ing jets at relativistic velocities. Both black holes as well as 
neutron stars have been identified as the compact, accret- 
ing object. The analogy of these systems to extragalactic 
quasars and Active Galactic Nuclei (AGN) make them ex- 
cellent laboratories for the study of the physics involved 
in accretion disks and ejection of relativistic jets associ- 
ated with accreting black holes (Mirabel et al. 19221 ■ The 
physics ruling black hole systems is essentially the same in 
galactic microquasars, hosting stellar black holes, and in 
AGN, with supermassive black holes, but the differences 
in time scales makes the study of some aspects in micro- 
quasars easier than in AGN. The characteristic timescale 
for the flow of matter onto a black hole is proportional to 
its mass, being of order of minutes in a microquasar of a 
few solar masses, but of thousands of years in a massive 
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black hole of 10 M Q . In addition, thanks to their prox- 
imity, microquasar jets have proper motions in the plane 
of sky with velocities about a thousand times faster than 
AGN, and two-sided jets can be observed. 

The microquasar GRO J1655-40 (X-ray nova Sco 1994, 
Zhang et al. [1994) was the second superluminal 
source discovered in our Galaxy, after GRS 1915+105 
(Mirabel & Rodriguez I1994[) . The two sources may also 
be peculiar in the sense that they show evidence suggest- 
ing that both systems contain a maximally spinning black 
hole (Zhang et al. ITM7|) . Radio images of GRO J1655- 
40 showed twin jets with apparent superluminal motion 
(Tingay et al. I1995f) moving in opposite directions at 
0.92c, and the distance was determined to be 3.2±0.2 kpc 
(Hjellming & Rupen 1995). Using the dust scattering halo 
observed by ROSAT, Greiner et al. I|1995[) determined a 
distance of 3 kpc, compatible with the previous determi- 
nation. Optical observations in 1996 provided an inclina- 
tion angle of 69.5±0.08°, a radial velocity semiamplitude 
of 228.2±2.2 km s -1 and a dynamical mass of the primary 
component of 7.02±0.22 M Q (Orosz & Bailvn ITW7|) . in- 
dicating that it is a black hole. Hubble Space Telescope 
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observations showed that the system moves in an eccentric 
orbit with a runaway space velocity (i.e., with respect to 
the Galactic rotation corresponding to its position in the 
Galactic plane) of 112±18 km s" 1 (Mirabel et al. l2UU2|) . 
which together with abundance anomalies found in the 
atmosphere of the donor star (Israelian et al. I1999f) pro- 
vide evidence for a supernova origin of the black hole in 
GRO J1655-40. 

ASCA observations of GRO J1655-40 in August 1994 
and August 1995 provided the first detection of absorption 
lines in an accretion powered source (Ueda et al. |I998). 
The energy of the lines was found to depend on the X- 
ray intensity, being 6.95 keV (Fe XXVI Ka) at 2.2 Crab, 
and 6.63 and 7.66 keV (Fe XXV Ka and K/3) at 0.27- 
0.57 Crab, revealing the presence of a highly ionized ab- 
sorber. Similar absorption features were also detected for 
GRS 1915+105 (Kotani et al. I2UUU|) . During the out- 
burst in 1997, ASCA observed GRO J1655-40 between 
25 and 28 February, when it was at 1.1 Crab, and re- 
vealed an absorption feature at 6.8 keV, interpreted as 
the blend of the two resonance Ka lines of Fe XXVI and 
Fe XXV (Yamaoka et al. I2UUT|) . Simultaneously, RXTE 
observations performed on 26 February 1997 showed emis- 
sion lines at 5.85 and 7.32 keV, interpreted as the first 
red- and blueshifted Fe Ka disk line in a Galactic source 
(Balucinska-Church & Church 2000), indicating the pres- 
ence of significantly ionized material in a region of the disk 
at ~ 10 Schwarzschild radii. 

The first observed X-ray outbursts of GRO J 1655- 
40 occurred in 1994-1995 (Zhang et al. IT5M1 
Harmon et al. 1995), and a 16-month-long outburst 
began on 25 April 1996 (Sobczak et al. [T533)l . After 
7 years of inactivity, GRO J1655-40 left quiescence 
again on 17 February 2005 (Markwardt & Swank |2"UU5|| . 
The X-ray evolution was followed with RXTE/ ASM 
(Homan et al. HzTTO and Swift (Brocksopp et al. EHI)- 
GRO J 1655-40 entered in February a low/hard state 
(Homan 2005) until it experienced a first outburst in 
March, moving into a high/soft state and reaching ~ 2 
Crab (Fig.^). The decay of the first outburst was followed 
by a month and a half of increasing X-ray flux and finally 
a strong outburst in very high state in May 2005, when 
the source reached more than 4 Crab. 

Here we present XMM- Newton observations performed 
in February-March 2005, aimed at obtaining detailed spec- 
troscopy during the rise to the high/soft state. We use the 
canonical distance of 3.2 kpc throughout (see note added 
in proof, however). 

2. Observations and data reduction 

We observed GRO J1655-40 with XMM-Newton some 
days before the rise to the soft/high state, on 27 February 
2005 (40 ks; obs. id. 0112921-301), and again on 14, 15 
and 16 March 2005, close to the maximum of the first 
outburst (15 ks each, obs. id. 0112921-401/-502/-601; see 
Fig. Hand Table^l. The EPIC-pn camera was used in the 
Burst mode, a fast read-out mode that permits the obser- 
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Fig. 1. RXTE/ASM light curve of GRO J1655-40, with 
dates of our XMM-Newton observations, and ASM hard- 
ness ratios: HRl, ratio of the ASM count rates (3- 
5 keV)/(1.5-3 keV), and HR2, (5-12 keV)/(3-5 keV). 
Data from the quick-look results provided by the 
ASM/RXTE team. 



vation of very bright sources with a high time resolution 
(Kuster et al. 1999). In the Burst mode, only one of the 12 
CCDs of the camera is used. The electrons released in the 
area of the CCD where the source is positioned are quickly 
shifted towards the read-out node to prevent pile-up, and 
then normally read out. In this mode, the spatial informa- 
tion in one dimension (shift direction) is lost and only 3% 
of the total time is effectively used, but bright sources can 
be observed without photon pile-up problems. The two 
MOS cameras were switched off for telemetry limitations, 
while the RGS spectrometers were operating in normal 
spectroscopy mode. 

XMM-Newton data were reduced using SAS 6.1, and 
XSPEC 11.3 was used for spectral analysis. Solar flares 
affected partially our observations, and some time inter- 
vals could not be used for analysis (in particular, the sec- 
ond half of the first observation lasting 20 ks, the first 2 
ks of the second and the first 5 ks of our last observa- 
tion). For the EPIC-pn data, only single events with en- 
ergy higher than 0.4 keV were taken into account. Events 
with RAWY>140 were excluded to avoid direct illumina- 
tion by the source, and the source spectra were extracted 
from a box including 10 RAWX columns at each side of the 
source position and ignoring also RAWY<10. The high 
luminosity of GRO J 1655-40 at the time of our XMM- 
Newton observations (almost twice as bright as the Crab, 
the source used for Burst mode calibration) made evi- 
dent a rate-dependence of the Charge Transfer Efficiency 
(CTE) not included in the calibration. We determined this 
dependence and corrected our spectra accordingly (see ap- 
pendix) . 

In the RGS spectra, the energy range 0.3-2.0 keV was 
initially included. We applied a preliminary version of 
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wavelength and time dependent corrections to the RGS ef- 
fective area, which will be implemented in the next release 
of the XMM-Newton software. Nevertheless, during our 
three March observations, strong differences were found 
between the responses of the different nodes of the RGS, 
probably due to the high flux of the source. We have thus 
ignored the energy ranges with evident problems, more 
precisely, energies above 1 keV for the first order in all 
cases, and between 0.85 and 1.2 keV in the second or- 
der of the RGS2 on the 14th March observation. Since 
both first and second order spectra of the two RGS are 
available, no energy range remains uncovered despite the 
rejection of these data. 

Simultaneous RXTE observations of GRO J1655-40 
are publicly available (Table ©. We have included the 
standard spectral products from the PCA and HEXTE 
instruments in our analysis, to constrain the model com- 
ponents at hard X-rays. When EPIC-pn and PCA spectra 
of March observations are attempted to be fit simultane- 
ously, large residuals appear for the PCA spectrum for en- 
ergies lower than ~ 6 keV, in the region of the instrumen- 
tal Xenon edges. Since the accuracy of the PCA response 
function is limited by the poor modeling of these edges 
(Jahoda et al. I2005[) . it is possible that the high flux of 
the source in these observations is revealing some inaccu- 
racy in the calibration, as in the case of the XMM-Newton 
instruments. We thus only consider the PCA spectrum for 
energies above 6 keV. 

Neither orbital variations, like eclipses or periodic 
modulations, nor dips are observed in the light curve, and 
the flux of the source remains constant within 30%. 

3. Spectral analysis 

3.1. February observation 

During the first observation, on 27 February 2005, be- 
fore the start of the first outburst, the source was in a 
low/hard state, with the XMM-Newton plus RXTE in- 
struments showing a spectrum dominated by an absorbed 
power law with photon index r=1.48±0.01 (uncorrelated 
90% error for 1 degree of freedom). This is in agreement 
with the results from RXTE observations starting only 
two days later, in the period 1-6 March 2005 (r=1.4- 
1.6, Swank & Markwardt I2005[> . and about 10% flatter 
than reported by Brocksopp et al. (20(23) from Swift ob- 
servations on 6 March 2005 (r=1.72±0.03). The fit with 
a power law leaves excess residuals at low energies both in 
the EPIC-pn and the RGS spectra, but the fit improves 
when the spectrum is modified by the scattering of dust 
in the line-of-sight. Greiner et al. l|1995j> found in ROSAT 
HRI observations in 1994 a halo around GRO J1655-40 
due to the scattering of the X-rays by interstellar dust. 
We modify thus our spectral model with the 'dust' model 
available in XSPEC, which assumes that the scattered flux 
appears as a uniform disk on the sky (as for a source con- 
stant in time) whose size has a 1/E dependence and whose 
total flux has a 1/E 2 dependence. 



GRO J1655-40, 14 Mar 2005 
XMM-Newton and RXTE 




1 10 100 

channel energy (keV) 



Fig. 2. XMM-Newton and RXTE observed spectra (upper 
panel), and unfolded spectra showing the contribution of 
the thermal disk emission (soft energies) and the power 
law (hard energies) to the total spectrum (lower panel). 
Residuals for the three March observations are shown in 
Figs. 12 Hand El 

The RXTE/PCA spectrum also shows residuals be- 
tween 6-7 keV, which are well fitted with an iron fluores- 
cence line with energy fixed at 6.4 keV, for which we derive 
an equivalent width of 135±20 eV. The line is not seen in 
the EPIC-pn data, and fixing the energy at 6.4 keV, we 
obtain an upper limit of 140 eV for the equivalent width. 
The presence of an iron fluorescence line in the low/hard 
state is confirmed by Swift observations one week after our 
XMM-Newton observation (Brocksopp et al. l2006[> . All pa- 
rameters of the fit for this first observation are listed in 
Tables H and 

3.2. March observations 

The three March observations have very similar spectra, 
with the XMM-Newton continuum spectra dominated by 
an absorbed multicolor disk (MCD), with kTi n increasing 
with time from 1.25 to 1.35 keV, and the higher energy 
spectra from RXTE represented by a soft power law, with 
r=2. 1-2.3 keV (see Fig.EJand Table©. The simultaneous 
fit of the EPIC-pn, RGS, PCA and HEXTE data with the 
MCD plus power law model leaves large residuals again 
below 0.6 keV in both the EPIC-pn and the RGS spectra, 
as well as in the range 8-15 keV in the XTE/PCA spectra. 

The soft excess in the XMM-Newton instruments, as 
in the case of the February observation, is most probably 
related to the dust scattering. Although the global fit im- 
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Table 1. Observation log 



Instrument 


Date 


Start 


End 


Exposure 


Energy range 


Flux (xlO ) 






(U lO) 


(U lO) 


(ks) 


[Ke V ) 


(erg s ^cm ^) 


XMM/EPIC-pn 


27-02-2005 


07:56 


19:34 


41.8 


0.4-10.0 


2.5 


XMM/RGS1 & 2, order 1 


27-02-2005 


07:46 


19:35 


42.5 


0.3-2.0 


0.5 


XMM/RGS1 & 2, order 2 


27-02-2005 


07:46 


19:35 


42.5 


0.54-2.0 


0.4 


RXTE/PCA 


27/28-02-2005 


20:37 


02:53 


13.0 


6.0-25.0 


5.9 


RXTE/HEXTE 


27/28-02-2005 


20:37 


02:17 


4.2 


20-200 


21 


XMM/EPIC-pn 


14-03-2005 


16:27 


20:36 


14.9 


0.4-10.0 


214 


XMM/RGS1 & 2, order 1 


14-03-2005 


16:17 


20:37 


15.6 


0.3-1.0 


2.9 


XMM/RGS1 & 2, order 2 


14-03-2005 


16:17 


20:37 


15.6 


0.48-1.8 


40 


RXTE/PCA 


14-03-2005 


13:02 


19:17 


12.9 


6.0-25.0 


54 


RXTE/HEXTE 


14-03-2005 


13:02 


19:17 


4.0 


20-200 


36 


XMM/EPIC-pn 


15-03-2005 


16:27 


20:36 


14.9 


0.4-10.0 


208 


XMM/RGS1 & 2, order 1 


15-03-2005 


16:17 


20:37 


15.6 


0.3-1.0 


2.9 


XMM/RGS1 & 2, order 2 


15-03-2005 


16:17 


20:37 


15.6 


0.48-1.8 


40 


RXTE/PCA 


15-03-2005 


14:08 


19:26 


11.6 


6.0-25.0 


48 


RX IE/HEX IE 


15-03-2005 


14:08 


19:36 


3.5 


20-200 


16 


XMM/EPIC-pn 


16-03-2005 


16:14 


20:23 


14.9 


0.4-10.0 


260 


XMM/RGS1 & 2, order 1 


16-03-2005 


16:04 


20:24 


15.6 


0.3-1.0 


3.2 


XMM/RGS1 & 2, order 2 


16-03-2005 


16:04 


20:24 


15.6 


0.54-1.8 


45 


RXTE/PCA 


16-03-2005 


17:08 


22:37 


11.1 


6.0-25.0 


70 


RXTE/HEXTE 


16-03-2005 


17:08 


22:37 


3.5 


20-200 


34 



proves after adding the dust halo effect, the RGS spectra 
still show an excess for energies lower than 0.5 keV. Using 
different abundance tables or leaving the abundances as 
free parameters did not help to fit this soft excess. This 
is probably related to the more limited sky area included 
in the RGS than in the EPIC-pn extraction regions. A 
detailed analysis of the halo, including spatially resolved 
spectroscopy, is currently underway and will be reported 
separately. As a first approximation, fixing the dust halo 
size to 10 times the beam size, from the fit with the EPIC- 
pn data we obtain in all cases a scattering fraction at 1 keV 
of 17-20%, similar to the 21% relative intensity of the halo 
found by Greiner et al. (^995). 



The quality of the fit in the 8-15 keV energy range 
spectra of the RXTE instruments during the three March 
observations improves with the addition of a smeared 
edge at 8.6-8.8 keV (see Table |3J|. An iron edge at 
8.81 keV was also detected in the 1994 outburst from 
ASCA observations (Ueda et al. I1998[> and at 8 keV 
with RXTE observations during the 1996-1997 outburst 
(Sobczak et al. I1999f) . Nevertheless, for the 15 and 16 
March observations, even after including the edge, the fit 
of the PCA spectra below 20 keV is poor. It is worth notic- 
ing that Sobczak et al. ilH)9!)i did also not find a good fit 
for the PCA spectra between 2.5 and 20 keV with MCD 
plus power law model, neither including emission lines, 
edges nor other effects, for the observations performed 
during the rise of the source to the first burst in 1996, 
in a phase of the light curve similar to our March obser- 
vations. 



3.3. Absorption features in the March observations 

Clear structured residuals appear also in the EPIC-pn 
spectrum, at ~ 6.8 and ~ 7.9 keV, corresponding to 
absorption features in the Fe K region (see Fig. [3] and 
Tabled}. In all three March observations the energy of the 
lines is the same, corresponding either to blended Fe XXV 
and Fe XXVI Ka and K(3 lines, or to blueshifted Fe XXV 
lines. In the latter case, given the energy resolution of the 
EPIC-pn (150 eV @ 6.4 keV), the velocity is not well con- 
strained. 

In the three March observations, the fit of the RGS 
spectra around the interstellar neutral oxygen edge at 
22.6-23.2 A (534-549 eV) with the phabs absorption 
model is unacceptable unless the "Vern" cross-sections 
(Verner et al. H996fl available in XSPEC are used. Even in 
this case, the fit leaves clear residuals, indicating the pres- 
ence of an additional absorption line at 23.5 A (Fig. ^J, 
which corresponds to the ls-2p interstellar absorption line 
of neutral oxygen. Given the energy resolution of the RGS, 
the observed maximum width of the line is compatible 
with it being a narrow line originating in the interstellar 
medium. The RGS spectra show also clear absorption fea- 
tures at 17.2 A and 17.5 A (Fig. 03), which correspond to Fe 
L2,3 (2pi/2, 2p 3 / 2 ) edges, also observed in the interstellar 
absorption towards some other sources (Wilms et al.|2000 
Paerels et al. 120011 Schulz et al. 12002(1 . As in these previ- 
ous cases, simple edge models do not provide a good fit 
for these features. 

Since both the line at 23.5 A and the edges at 17.2 A 
and 17.5 A are caused by interstellar gas but are not well 
reproduced by the simple phabs model in XSPEC, the 
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Parameters 27 Feb 2005 14 March 2005 15 March 2005 16 March 2005 



N H (xl0 21 ) cm" 2 




5.2±0.2 


5.46±0.02 


5.57±0.02 


5.45±0.03 


Dust fraction 51 




0.15±0.03 


0.186±0.002 


0.169±0.003 


0.184±0.003 


kTin (keV) 






1.231±0.003 


1.268±0.003 


1.326±0.003 


(km) 






14.25±0.05 


13.25±0.05 


13.55±0.05 


r 




1.48±0.01 


2.10±0.01 


2.27±0.01 


2.13±0.01 


K c (ph keV" 1 cm- 2 s- 1 @ 
L d o.4--io kcv (ergs 


1 keV) 


0.0385±0.0007 


0.94±0.02 


0.76±0.07 


1.05±0.02 


3.7xl0 35 


3.4xl0 37 


3.4xl0 37 


4.2 xlO 37 



a Dust scattering fraction at 1 keV. Halo size is fixed to 10 times the detector beam-size at 1 keV. 

b For a distance of 3.2 kpc and an inclination of 70°. 

c Normalization constant for the power law component. 

d Luminosity in the range 0.4-10 keV, for a distance of 3.2 kpc. 



Table 3. Features in RXTE/PCA spectra 



Parameters 



27 February 2005 14 March 2005 15 March 2005 16 March 2005 



Line energy (keV) (FIXED) 
la line width (eV) 
Equivalent width (eV) 
Identification 

Edge energy (keV) 
Edge optical depth 
Edge width 



6.4 
300±100 
135±20 
Fe I 



8.68±0.05 
0.44±0.05 
1.2±0.2 



8.57±0.03 
0.5±0.2 
1.9±0.5 



8.77±0.04 
0.57±0.05 
1.3±0.2 



recent new version of the TBabs model (Wilms et al. 2006 
http://astro.uni-tuebingen.de/~wilms/research/tbabs/) 
has been used. This has been shown to provide a good fit 
to the complex oxygen edge (Juett et al. l2004[) . including 
two edges and five absorption lines, and to the iron absorp- 
tion features (Juett et al.|2006). It provides a good fit also 
for our RGS features (see Figs. 0] and [5]), with a hydrogen 
column N H : 7.3(±0.3), 7.0(±0.3), 6.9(±0.5) x 10 21 cm- 2 , 
and an overabundance of oxygen with respect to solar of 
1.45(±0.15),1.6(±0.1),1.45(±0.15) (for days 14, 15 and 
16 March 2005 respectively). 

4. Discussion 

4.1. The dust halo 

The effect of the scattering dust halo is clearly detected in 
the soft X-ray spectra. The scattering fraction at 1 keV is 
changing during the three March observations, decreasing 
from 18.6(±0.2)% to 16.9(±0.3)% in the first 24 hours, 
and then increasing again to 18.4(±0.3)% during the next 
day. This temporal variation is consistent with a constant 
dust halo but a variable X-ray source flux: the time needed 
for the soft photons originating on the disk to travel to the 
scattering site far away from the source causes a delay in 
the flux variations of the dust halo with respect to the vari- 
ations of the source. As a result, just after an increase of 
the disk temperature and thus of the source soft spectrum, 
the fraction of scattered photons with respect to the source 



is temporally smaller (during the second March observa- 
tion) until the increasing source photons reach the outer 
part of the halo and their scattered photons also reach the 
observer, increasing again the scattering fraction (third 
March observation). This effect has been known for some 
time. It was also observed by ROSAT and used to estimate 
the distance to GRO J1655-40 (Greiner et al. ITMijl . 

4.2. Continuum 

The maximum disk temperature observed in the March 
observations, increasing from 1.25 to 1.35 keV, is con- 
sistent with the values commonly observed in black 
hole binaries in high/soft state (see for instance 
McClintock & Remillard 2006), and with the Swift obser- 
vations of the same period (Brocksopp et al. I2006fl . The 
inner radius, R; n ~ 13 — 14 [D/3.2kpc] km, is only a lit- 
tle larger than the gravitational radius of the black hole, 
R g = GM/c 2 = 10 km, consistent with a rapidly spinning 
black hole (Zhang et al. HW7jl . 

From the fit parameters of the MCD model the lumi- 
nosity of the disk can be calculated as L x = 47rR 2 n oT- t n 
(Makishima et al. I1986|) . With our values, the accretion 
disk luminosity is 6 — 7 x 10 37 [D/3.2kpc] erg s _1 . 

With the accretion disk luminosity, the ac- 
cretion rate can be calculated as 2L x R; n /GMg 2 
(Makishima et al. [T^grJjl . where g = (1 - Rg/Rn) 1 / 2 
is a correction for the general relativity, R g is the 
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Table 4. Absorption features in EPIC-pn spectra 



Parameters 


14 March 2005 


15 March 2005 


16 March 2005 


Line energy (keV) 


6.82±0.06 


6.71±0.06 


6.75±0.05 


la line width (eV) 


noljg 


200±100 


250±50 


Equivalent width (eV) 


60lf 


100±60 


110±30 


Identification 51 


Fe XXV Ka (E 6.700 keV)/Fe XXVI Ka (E 6.952/6.973 keV) 


v (km s ) 


-5200±2600 


-400±2700 


-2200±2200 


Line energy (keV) 


8.04±0.06 


7.85±0.15 


7.9±0.1 


lcr line width (eV) 


lOOir 


700±200 


600±200 


Equivalent width (eV) 






240±40 


Identification 51 


Fe XXV K/3 (E 


7.872/7.881 keV)/Fe 


XXVI K/3 (E 8.246/8.252 keV) 


v (km s" 1 ) 15 


-5900±2200 


1200±5700 


-700±3800 



a Based on CHIANTI data base (Dere et al. \TM7\ Landi et al. 1211061 
b Blueshift with respect to Fe XXV lines. 
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Fig. 3. EPIC pn residuals of GRO J1655-40 spectra of 
our three March 2005 observations at zero velocity, after 
fitting an absorbed multi-temperature disk model. 



gravitational radius and M is the mass of the black 
hole. From the accretion disk luminosity, the disk inner 
radius and the black hole mass of GRO J1655-40, we 
obtain an accretion rate during the high/soft state of 
- 10~ 8 M yr" 1 . 

The rise of the disk blackbody component during the 
high/soft state in March is reflected in the power law 
index of the corona, which increases from 1.48±0.01 in 
February 27 to values higher than 2 in the March obser- 
vations. With the rise of the thermal disk emission, the 
soft luminosity increases and the electrons in the corona 
are cooled more efficiently, resulting in a softer spectrum 
of the comptonizcd photons, i.e., a steeper power law. 



4.3. Interstellar absorption 

The interstellar OI Ka line is clearly detected at 23.5 A 
in our three March 2005 RGS spectra. This line was 
also found in the RGS spectra of other sources, and 
distinguished from the instrumental components around 
the interstellar oxygen edge, at 23.05 and 23.35 A (de 
Vries et al. EMiSj) . Juett et al. (2004) detected the same 
line in the high resolution spectra of seven X-ray binaries 
using the Chandra/HETGS, as part of their study of the 
structure of the oxygen absorption edge caused by the in- 
terstellar medium. More recently, it has also been detected 
in the Cyg X-2 spectrum (Costantini et al. I2005J) . in the 
observations that have allowed the first spatially resolved 
spectroscopic study of a scattering dust halo, as well as in 
the absorption towards LMC X-3 (Wang et al. 120051 and 
the recently discovered black hole candidate XTE J1817- 
330 (Sala & GreinerEDDEJ- 

To check the interstellar origin of the OI Ka line ob- 
served in the RGS spectra of GRO J1655-40, we have 
looked for a correlation of the equivalent width of the line 
and the hydrogen column density, using the above men- 
tioned works (Fig. E| . We find that the equivalent width 
of the OI Ka line, EW, is indeed correlated with the col- 
umn density, N H , as EW(eV)=0.5(±0.1) + 2.8(±0.3) x 
10 22 N H (cm- 2 ). The values observed for GRO J1655-40 
fall well within this correlation for the ISM. We find that 
an alternative correlation with zero regression constant 
could be EW(eV)=3.6(±1.6) x 10 22 N H (cm- 2 ), but this 
would underestimate the equivalent width for the sources 
with hydrogen column density less than 5 x 10 21 cm~ 2 . 

4.4. Orbital dependence and inclination of the binary 
system 

Since the orbital period of GRO J 1655-40 is ~ 2.6 days 
and our three March observations were taken with one 
day intervals, we are covering approximately one whole or- 
bital cycle. We can extrapolate the Orosz & Bailyn ( 1997) 
ephemeris, to obtain the orbital phase of our observations 
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Juett et al. (2004): several sources 
de Vries et al. (2003): several sources 
Costantini et al. (2005): Cyg X-2 
Sala & Greiner (2006): XTE J1817 -330 
This work: GRO J1655-40 



5x10" 



N„(cm- 2 ) 



10- 



2x10* 



ivelength (Angstrom) 



Fig. 4. RGS1 (black) and RGS2 (grey, contributing only 
below 20A and above 24A.) around the oxygen edge in 
the first order spectra of the 14th March 2005 observation, 
fitted with the new TBabs model. 



Fig. 6. OI K-a line equivalent width versus hydro- 
gen column density for several galactic sources. The 
dashed line indicates an approximate linear correlation, 
EW(eV)=0.5 + 2.8 x 10- 22 N H (cnr 2 ). 




16.5 17 17.5 18 18.5 19 

channel wavelength (Angstrom) 



Fig. 5. First and second order (with lower count rate) 
spectra of GRO J1655-40 in the range 16-19A, obtained 
with RGS1 (black) and RGS2 (grey) on 14 March 2005, 
showing the ISM Fe L edges fitted with the new TBabs 
model. 

(error is 0.08 in all cases): 0.31-0.40 (27 February), 0.17- 
0.24 (14 March), 0.55-0.62 (15 March) and 0.93-0.00 (16 
March). 

Given the inclination of the system (between ~ 70°, 
Orosz & BailynEHnil van der Hooft et al. H33SI and - 85°, 
Hjellming & Rupen 1995), an orbital phase close to zero 
corresponds to having the donor star situated closer to the 
observer, while in phase 0.5 the disk would be in front of 
the secondary star. This means that if any of the observed 
spectral features were arising from the illuminated face of 
the secondary star, it would have its maximum at phase 
0.5 and would not be present at phase zero. This is not the 



case in our observations. Another possible orbital effect 
would be that the disk emission were absorbed by the 
stellar wind of the secondary star, which would produce 
increased absorptions close to phase zero, i.e. on 16 March, 
which is also not the case. 

Finally, neither dips nor eclipses were observed in the 
disk thermal emission at any of the phases. Given the 
size of the donor star, ~ 5R© (Orosz & Bailyn 1997), 
the binary separation (1.17 x 10 12 cm), and assuming that 
the soft X-ray emission originates in the central 200 000 
km of the disk (see below, section I4.5|l , the duration of a 
possible eclipse of the soft X-ray emission by the donor 
star would last more than 2.5 hours (i.e., a change in the 
orbital phase of 0.04) close to phase zero. This should have 
been clearly visible during our 16 March observation. 

The fact that no orbital modulation is observed in the 
X-ray spectra, provides an upper limit for the inclination 
of the system. With the parameters mentioned above for 
the sizes of the system, the inclination must be smaller 
than 73° (for the innermost 200 000 km of the disk sur- 
face to be visible in all orbital phases). This limit is in 
agreement with the inclination determined from the op- 
tical light-curve by Orosz & Bailyn <fT§9T|) . 69?5 ± 0?8, 
and by van der Hooft et al. lfT5Mj> . 67? 2 ± 3? 5, but in- 
compatible with the inclination inferred from the radio 
jets, ~ 84° (Hjellming & Rupen 1995). This points, as 
suggested by Orosz & Bailyn (1997), to an inclination of 
the jet axis of about 15° with respect to the normal to 
the orbital plane. Alternatively, non-symmetric jet ejec- 
tions could have lead to a wrong inclination inference 
(Hjellming & Rupen DM^ . 
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4.5. The Fe XXV/Fe XXVI absorber 

Clear highly ionized Fe absorption lines are detected in the 
EPIC-pn spectra of the three March observations in the 
high/soft state. It is worth noting that Ueda et al. I|1998[) 
found simultaneously Fe XXV absorption features and an 
iron edge, while they detected no edge when absorption in 
the iron K band corresponded to Fe XXVI. 

Assuming that the absorbing plasma is photoionized, 
the ionization state of the observed elements provides in- 
formation on the conditions and location of the absorber. 
The presence of He-like Fe ions indicates an ionization pa- 
rameter £ =L/nr 2 ~ 10 3 erg cm s _1 (Kallman et al. 1996). 
The lack of lower ionization absorption implies that the 
absorber cannot have a large extent. 

The equivalent width of the Fe lines is increasing be- 
tween the first and the last March observations. This may 
be pointing out either an increase in the total column den- 
sity, or in the ionization state of the gas. 

Ueda et al. l|1998fl determined the iron column den- 
sity of the plasma to be 10 19 — 10 20 cm' 2 from the 
observed equivalent width and the curve of growth of 
the Fe XXV Ka line, which relates the expected equiv- 
alent width to the iron column density. Using their curve 
of growth and assuming the detected features are only 
Fe XXV, our equivalent width of the Fe XXV Ka line (be- 
tween 50 and 160 eV) corresponds to a column density of 
10 19 and 5 x 10 20 cm -2 , which is similar to the values found 
by Ueda et al. <|1998fl . Assuming cosmic abundances, this 
corresponds to a hydrogen column density in the range 
2 x 10 23 — 10 25 cm -2 . Since the observed absorption lines 
may be a blend of Fe XXV and Fe XXVI, these must be 
taken as upper limits. In addition, for a hydrogen column 
density larger than 10 24 cm -2 , the line absorber would be 
optically thick to Thomson scattering (Ueda et al. 1998). 
We consider this to be unlikely and, taking into account 
the errors in the equivalent widths, the column density 
could be then in the range 2 x 10 23 — 10 24 cm~ 2 in all 
three observations. 

With the photoionization parameter, £ ~ 10 3 cm s , 
and the flux detected above ~ 9 keV (X-rays photoionizing 
He-like iron ions), L>gk v ~ 5 x 10 36 erg s" 1 (for a distance 
of 3.2 kpc), the column density indicates a distance to the 
central source between 50 000 and 200000 km. Assuming 
that the disk radius is 70% of the Roche lobe radius, i.e., 
4 x 10 6 km, the Fe XXV absorber extends to less than 5% 
of the disk surface. 

In the case that the features correspond to blueshifted 
Fe XXV lines, the error in the energy determination is too 
large as to reach any strong conclusion about its blueshift. 
Nevertheless, assuming the wind velocity is constant, the 
common range of blueshifts of the three March observa- 
tions is reduced to 2600-4500 km s _1 . From the parame- 
ters derived above for the Fe XXV absorber, and assuming 
constant density and spherical symmetry for the expand- 
ing wind, the mass loss rate 47rr 2 /jvb would be in the range 
(2 - 13) x 10~ 7 b M yr _1 , with b being the filling factor. 
This upper limit is a factor 20-130 larger than the accre- 



tion rate, which would indicate a highly non-stationary 
situation. However, there is considerable uncertainty in 
the wind mass loss, since both the wind density and wind 
velocity are likely a function of radius and the wind may 
be conical rather than spherically symmetric. Exploring 
these effects goes beyond the scope of this paper. 



5. Summary 

GRO J 1655-40 was observed with XMM-Newton during 
the low/hard state in February 2005 and in the soft/high 
state in March 2005. The continuum spectrum is domi- 
nated by an absorbed power law in the low/hard state 
and by a multicolor disk blackbody in the soft/high state, 
with the maximum temperature increasing during the 
high state. In addition, the X-ray spectrum during the 
March observations shows several spectral features: 

- Interstellar neutral oxygen causes an OI Ka absorp- 
tion line at 23.5 A. A correlation between the equiva- 
lent width of this line and the hydrogen column density 
has been obtained from published observations of several 
galactic X-ray sources, and the comparison of our observa- 
tions with that correlation confirms the interstellar origin 
of the OI Ka absorption line. 

- Fe XXV, possibly blended with Fe XXVI, is detected 
in absorption in the EPIC-pn spectra, indicating the pres- 
ence of a highly ionized absorber, with an ionization pa- 
rameter 1000 erg cm s , extending less than 200 000 km 
from the central source, corresponding to a 5% of the ac- 
cretion disk. 

Photoionized absorbers are seen in AGN as well, also 
with expansion velocities. The detection of photoionized 
winds in microquasars is then another similarity between 
these two kinds of systems. Nevertheless, winds in AGN 
are not uniform, but rather distributed in clouds of mate- 
rial. If this were the case in GRO J1655-40, the absorption 
systems would go in and out of the line of sight and cause 
temporal variations in the absorption features, which is 
not observed in our XMM-Newton observations. 
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Note added in proof: 

Within a few days of placing the refereed accepted ver- 
sion of this paper on astro-ph, two other papers appeared 
about the same source GRO J1655-40, which partially af- 
fect our conclusions and are commented on below: 

1. Foellmi et al. l|2tT06l have recently obtained a new de- 
termination for the distance to GRO J1655-40 and find 
an upper limit of 1.7 kpc. In this case, the luminosity 
in the range 0.4-10 keV listed in Table [21 would be 
1.0 x 10 35 , 9.6 x 10 36 , 9.6 x 10 36 and 1.2 x 10 37 ergs" 1 
(for days 27 Feb, 14, 15 and 16 March respectively), 
and the disk inner radius would be only 7.6 km, 7.0 km 
and 7.2 km (for days 14, 15 and 16 March 2005). We 
note, however, that these values are smaller than the 
gravitational radius of the black hole, R g = GM/c 2 ~ 
10 km. Therefore, either the distance to GRO J1655- 
40 is at least 2.5 kpc, or the model for deriving the 
inner radius of the accretion disk requires improved 
modifications. 

2. Miller et al. I|2006fl reported on a Chandra observa- 
tion of GRO J1655-40 performed on 1 April 2005, two 
weeks after our last XMM-Newton observation. The 
HETGS spectrum shows 90 absorption lines signifi- 
cant at the 5cr level of confidence or higher. From all 
of these, only the Fe XXV Ka and Kf3 lines are also 
present in our March observations. We have computed 
the upper limit for the flux of some of the strongest 
lines detected by Chandra in the RGS energy range, 
in particular, for the Fe XXIV 5p, 4p and 3p lines 
at 7.159, 7.979, 10.604, 10.649 A. Fixing the energy 
and FWHM of the lines to the values obtained by 
Miller et al. ( 2006 ) , we obtain the following upper lim- 
its at 3<7 confidence level for the flux of the lines: 4.2, 
0.36, 1.4 and 1.3(xl0 _3 phcm _2 s _1 ), for the observa- 
tion on 15 March 2005. These upper limits are between 
a factor 3 and 30 smaller than the fluxes measured in 
the Chandra spectrum and suggests a change in the 
ionization state between the two observation periods. 

Appendix A: EPIC-pn CTE correction 

An inaccurate calibration of the EPIC-pn charge transfer 
efficiency (CTE) leads to a slightly wrong energy deter- 
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Fig. A.l. Calibration of the gain linear factor / as a 
function of rate per pixel r. 

ruination of a detected photon, which becomes evident in 
the large residuals around sharp features like the instru- 
mental Si and Au edges. From our GRO J1655-40 ob- 
servations, the shift in energy has been found to be rate 
dependent, being stronger at the center of the PSF, which 
implies that it can not be directly corrected in the ex- 
tracted spectrum. We determined the improvement of the 
CTE (formally handled as a linear gain factor) for differ- 
ent rates, selecting and evaluating the energy gain linear 
factor for spectra extracted from different regions of the 
detector (Fig lA.lj) . Crab calibration observations in the 
Burst mode were also included to improve the determina- 
tion of the dependence of the gain linear factor / with the 
rate per pixel r (cts/s/pixel), which we find can be approx- 
imated by f = 0.98 + 0.015r- 2.2 x 10 _3 r 2 + l.l x 10~ 4 r 3 . 

We used this linear gain to correct the energy of our 
event files: for each photon, the PI channel was divided by 
the gain factor / according to its rate r, assigned to each 
photon depending on the RAWX column. To avoid the 
discretization resulting from the pixelization of the PSF in 
the RAWX columns, the r values were randomized within 
each RAWX column. 

After correcting the event tables with this linear gain, 
no more residuals appear around the Au and Si edges in 
the new extracted spectra. We confirmed that the correc- 
tion is also valid at the high energy range, by checking 
that, after the correction, the spectral features at 7-8 keV 
have the same central energy both at the center of the 
PSF (at high rate) and at the wings (low rate). 




